White-light parametric instabilities in plasmas 
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^ — ' Parametric instabilities driven by partially coherent radiation in plasmas are described by a gen- 

, eralized statistical Wigner-Moyal set of equations, formally equivalent to the full wave equation, 

■ coupled to the plasma fluid equations. A generalized dispersion relation for Stimulated Raman 

£\j ' Scattering driven by a partially coherent pump field is derived, revealing a growth rate dependence, 

with the coherence width a of the radiation field, scaling with 1/a for backscattering (three- wave 
O i 1 process), and with 1/a 1 / 2 for direct forward scattering (four-wave process). Our results demonstrate 

the possibility to control the growth rates of these instabilities by properly using broadband pump 
radiation fields. 
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Parametric instabilities are pervasive in many fields of science, associated with the onset of nonlinear and collective 
effects such as solitons, vortices, self-organization, and spontaneous ordering. Recent developments in light sources 
and laser technology continue to reveal novel features of the parametric instabilities, for instance in nonlinear optics, 
with the recent experimental discovery of white light solitons or in plasma physics, in the realm of relativistic 
nonlinear optics ■ The standard theoretical approach to study parametric instabilities is based on a coherent wave 
description which is clearly limited because, in most systems, waves are only partially coherent, with incoherence 
either inherently induced by fluctuations, or induced by external passive systems (e.g. random phase plates in inertial 
Oh! confinement fusion (ICF)). Recent theoretical work in nonlinear optics, triggered by the work of Segev and co-workers 
C/3 ' [H, led to the development of techniques capable of describing the propagation and the modulation instability of 
partially coherent /incoherent "white" light in nonlinear media Q. The critical underlying assumption of all these 
C/3 ■ models is the paraxial wave approximation, valid in transparent media for radiation beams not tightly focused, which 
reduces the problem of electromagnetic wave propagation in dispersive (and diffractive) nonlinear media to the search 
of a forward propagating solution, formally described by the nonlinear Schrodinger equation. While in nonlinear 
optics, and for the conditions studied so far, such approximation is clearly valid, in plasma physics it is not. The 
instabilities associated with the partially reflected backscattered radiation Q are critical in many laser-plasma and 
astrophysical scenarios Q , and the paraxial approximation has limited applicability even in the description of forward 
I ■ scattering instabilities driven by ultra intense lasers in underdense/transparent plasmas 

Inclusion of bandwidth/incoherence effects in laser driven parametric instabilities in plasmas and in three- wave 
processes is a long-standing problem 0, Q , because incoherent pumps can decrease the growth of the laser driven 
instabilities in ICF, Fast Ignition, or novel laser amplification schemes The difficulty resides in the lack of an 
■<^J- ■ appropriate theoretical framework where a statistical description of the radiation is natural. The Wigner-Moyal 
(WM) formalism of quantum mechanics [l(| provides a natural path to build such a statistical description of the 
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radiation 11]. However, it is important to point out that the standard WM approach is valid only for Schrodinger 
like systems, where backscattered radiation is neglected. Previous attempts 12J have only been able to describe direct 

^ • forward stimulated scattering, failing to describe in a general way parametric processes driven by "white" light in 
plasmas. To overcome this difficulty, we have recently developed a generalized WM statistical theory of radiation 
1 1 31 ] . or generalized photon kinetics (GPK), formally equivalent to the full wave equation, valid for partially coherent 

C3 ' electromagnetic wave propagation in nonlinear dispersive and diffractive media. 

In this Letter, GPK is employed to derive the general dispersion relation for Stimulated Raman Scattering (SRS) 
driven by a spatially stationary radiation field with arbitrary statistics, thus valid for all ranges of coherence of the 
pump field i.e. from a coherent plane wave pump to an incoherent pump. In this dispersion relation both three- wave 
processes and four-wave processes are considered, and for a plane wave pump field the standard results are recovered 
@, EH- Analytical results are derived for different regimes of SRS and wavenumber ranges, showing universal 
decays of the growth rate with the bandwidth er, with a 1/a dependence for backscattering, and a slower decay, with 
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l/cr 1 / 2 , for forward scattering. Numerical solutions confirm the theoretical predictions, highlighting some of the most 
important consequences of white light in SRS. 

In our discussion, we use normalized units, such that length is normalized to c/u> p0 , time to 1/lu p q, mass and 
charge to the electron mass m e and the electron charge e, respectively, and where c is the velocity of light in vacuum, 

1 /2 

and ojpQ — (4we 2 n e o/m e c 2 ) . We model the plasma as a cold uniform electron fluid with a fixed ion density 
^io = n e0 (in normalized units rt^o — n e0 = !)• The normalized vector potential of the circularly polarized pump field 
a p — eA p /m e c 2 is described by a p (r, t) — 2~ 1 l 2 {z + iy)a$ j dkA(k) exp [i(k ■ r — w(k)t)], where w(k) is the dispersion 

1/2 

relation for plane circularly polarized monochromatic waves in a uniform plasma, oj(k) = (k 2 + I/70) , and ^4(k) 
can include a stochastic phase dependence ip(r, t), as A(k) = A(k) exp (iip(r, t)). The only restriction on the form of 
a p (r, t) is that the Klimontovich statistical average of the two point correlation function, (a*(r+y/2, t)-& p (r— y/2, i)) = 
a§m(y), is independent of r with m(0) = 1, and |m(y)| is bounded between and 1 i.e. the field is spatially stationary. 
This restriction is introduced only because of the perturbation technique we are employing; the formalism described 
here is valid for any field dependence. When the uniform plasma is irradiated by the electromagnetic field, described 
by the normalized vector potential a = a p + a, the normalized electron density n e = 1 + n satisfies [l5| : 

# + ^)ft=^V?(<Re|Va])) (1) 



with 70 = y 1 + ( a p • a p) = \J 1 + a o ' wriere (• ••) denotes a statistical average, and ~ denotes first-order quantities. 
Usually, the driving term on the right hand side of Eq. |T|) is described with the standard approach based on the 
wave equation for the vector potential 0, Il5| . However, this technique does not allow for the study of "white" light 
parametric instabilities. GPK can address the general two mode problem. The radiation field a is described by 
two fields <j>,X — (a ± i^/7o<9ta)/2, in terms of which the wave equation can be written as two coupled Schrodinger 
equations [13(. Introducing four real phase-space densities Wo = W^— W xx , W\ = 21m [W^, x ], W2 = 2Re [W^, x ], W3 = 
+ W xx , where the Wigner transform Wf. s is defined as Wf. g (k, r, t) = (^r) 3 / e* k y f* (r + ^, i) ■ g (r — |, t) dy, 
it is possible to derive a set of transport equations for Wi [13|: 

d t W +£(W 2 +W 3 ) =0 (2) 

d t W 1 - g{w 2 + w 3 ) - w 2 = (3) 

\/7o 



d t W 2 - CW + QW X + Wi = (4) 



d t W 3 +CW -GWx = (5) 



where the operators C and Q obey: 



C ■■ - x k ■ V r x ( ^ sinQv r ■ V k j ((>) 




7b \*--r +V7^ - cos h;Vr-V k (7) 



with the arrows denoting the direction of the operator, and sin(...) and cos(...) represent the equivalent series expansion 
of the operators; we observe that the left arrow operator V r acts on ^/to (^J , while the right arrow operators (Vk> 

V rj V r ) act on W%. Equations ([2E]) are formally equivalent to the full wave equation for a in a plasma. 

In order to close the system of Eqns. ([1] [MS]), and to determine the corresponding dispersion relation, it is nec- 
essary to linearize Eqns. (J2HS1) , noting that up to first order Wq = y/jop (k) w(k) + Wo(k, r, t), W\ — Wi(k,r,t), 

W 2 = -p (k)^f- + W 2 (k, r, t), W 3 = po(k) (l + 3S^) + M> 3 (k, r, t), where p (k) = W ap . ap is the zero-order photon 
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distribution function. In analogy with the standard techniques in plasma physics, po(k) can be thought of as the 
equilibrium distribution function of the photons. Furthermore, we observe that W2 + W3 = 2WR e [ a .5]. Linearization 
of Eqns. ([2]-[5]), followed by time and space Fourier transforms (d t — > uil, V r — * —ik L ), leads to: 



Rc [a„-al 



= 2 T 




(8) 



where £)^(k) = lj 2 l =F 2 [k • — cj^ lj (k =p 2 )] ' an< ^ w ith •T'lfl'W.kL denoting the Fourier transform of <?(r, t). In 
order to obtain the plasma response ( ~Y the same technique is followed, leading to 
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1 LF[Re [a p - a]; 



(9) 



Integrating Eq. (|8|) in k, and using J Wf. s dk = f * ■ g, Eqns. (|8|9p can then be combined to give the exact dispersion 
relation for electron plasma waves in the presence of broadband radiation: 
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with D ± (k) — (u>(k) iwi) 2 — (kik^) 2 — Equation (fTU)l is the central result of this paper, and it generalizes 
the seminal result of Decker et al [l5| for pump fields with arbitrary statistics. It can be interpreted as the statistical 
average of £>+^ + D-(k) over the distribution of photons. For a pump plane wave, with wavenumber ko, po(k) = 
a 2 (5(k — ko), and Eq. (fTOf leads to the same dispersion relation as derived in Ref. [HI]. Recently, a dispersion 
relation with two pump waves was also obtained [1 61 ] , which also can be derived from Eq. (|10p for two photon beams 
p (k) =a^ 1 5(k-koi) + ag 2 5(k-k 02 ). 

To illustrate some of the most important consequences of white light in SRS, we consider the one dimensional 
scenario, for a water-bag zero-order photon distribution function powb(^) = Oo/( cr i + &2)(®(k — ^0 + f 1) — — 
ko — 02)), where 0(fc) is the Heaviside function. With this choice for po, several analytical results can be derived, 
highlighting the influence of white light in parametric instabilities. For powb(^), the random phase ip(x) is such 
that the autocorrelation function satisfies (exp (— iip (x + % ) + tip (x — §))) = exp (— iyE) sin (ya)/(ya), with W — 
((T2 + ci)/2 and k = (ko + (02 — cti)/2). The correlation length of this distribution is ~ 7r /V2a. The dispersion 
relation (I10|) for this distribution function, valid for all values of ko, ao and 171,2, is: 
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where 
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2kl(uL + kL)y/Qo {2a + ujqi — CJ02) / (Qokl — Q ± (wl + fc^) 2 ). It is now possible to study the effect of a broadband 
photon distribution on stimulated Raman forward scattering (RFS), relativistic modulation instability (RMI) or 
stimulated Raman back scattering (RBS). 

Analytical results can be obtained in the case of an underdense medium 1 /70 <C ko — ai, which also guarantees that 
ko > The first condition states that the medium is underdense for all the photons in the distribution, while the 
second assures that Powb(^) represents a broadband source of forward propagating photons. We further observe that 
under this approximation no order relation for 02 needs to be assumed. The underdense approximation is equivalent 
to neglecting the arctanh terms in (jlip . since cuoi ~ fco + (— l) l 0ij and w 0. We observe that by neglecting I/70 
when compared with ko, the simplified dispersion relation is still valid for RBS, but we loose the ability to capture 
RMI. For RFS, we then Taylor expand the log term in Eq. (|11[) . in the underdense approximation feo — cr± 1/^/tq. 
The resulting polynomial equation can the be evaluated near the wavenumber k^ for maximum growth rate, such 
that kff RFS s» 1/^/70, and lul ~ l/v/7b + S, with 6 -C 1, to yield the maximum growth rate for RFS, Trfs = Im[5]: 



rRFS = 



2\/27 V( fc o-^i) (ko + a 2 ) 



(12) 
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In the limit of — > 0, the standard monochromatic result, valid for all intensities, is obtained [l5|], with Trfs °c o-q 1 
for ao 3> 1. The effect of the bandwidth on four wave processes, such as RFS, predicted by Eq. (p~2|) is qualitatively 
different from the effect of the radiation bandwidth in 3-wave processes, such as RBS. Eq. (|12|) shows that Trfs 
increases(decreases) for increasing <7i(cx2); this can be interpreted from the monochromatic result as due to the 
decrease(increase) of the average wavenumber of the distribution of photons. For RBS, the same technique can be 
followed, but now must be resonant (corresponding to the contribution of the downshifted photons of the highest 
wavenumber photons in the distribution function), with ujl ~ 1 /s/^Yo "F *Frbs > with Trbs <C 1/\/7o as in the usual 
treatment of RBS [3], away from the strongly coupled regime, which means that the instability occurs in regions close 
to kff RBS ps 2(fco + a 2) — Furthermore, Df ~ 0, corresponding to the contribution of the downshifted photons 

of the lowest wavenumber photons in the distribution function, establishes the lower limit of the range of unstable 
wavenumbers given by fc™RBS ~ 2(fco — (X\) — 1/^/to. The maximum growth rate for RBS is: 

iral k a + g 2 1 

1rbs = — rjK : —2 — 

870 1 "f" 5°/ 2 ? 



where we have also assumed in the derivation of Eq. (|13[) that Trbs < 01 + 02- In the opposing limit, as <Ti 2 - * for the 
plane wave limit, we obtain Trbspw = aoV^o/v^To^ 4 , thus recovering the standard result 0]. For ao S> 1, Eq. (|13p 
shows the scaling Trbs oc a ^ 4 as in the monochromatic case (T^]. A comparison between Eq. (fT2j) and Eq. (fl3|) shows 
the stronger dependence of RBS on the bandwidth of the radiation, for cr 2 < fee In fact, for fixed fco, ao, and ai, RFS 
scales with cx 1/^/02, while RBS goes as cx 1/<J2- We have found this behavior for other distribution functions (e.g. 
asymmetric lorentzian/gaussian distribution of photons), and, within the aproximations discussed here, it is possible 
to show that the scaling of RBS and RFS with 02 is independent of the exact shape of the distribution function. 
The wavenumber for maximum growth kff is independent of 172 for RFS, while for RBS it depends linearly on 02. 
Furthermore, Eq. (fl~3|) also shows that for large values of 02, the growth rates for RBS satures at Fp a Bg = iraQ/8j^ 2 . 

To illustrate these features, we have solved Eq. (fTTj) numerically, for the range of fc_L required to capture RFS, RBS, 
and RMI, and for different 02, keeping fco, ao, and 01 constant. In Figure 1, the dependence of the maximum growth 
rate with 02 are shown, demonstrating that RFS is less sensitive to 02 than RBS, as predicted by the theory. The 
theoretical curves, given by Eqns. (|12ll3p , can not be distinguished from the numerical solution. RMI depends more 
strongly on 2 than RFS, but still decreases more slowly than RBS, with Prmi/F(02 = 0) ~ 1/ (l + 1.35(02/fco) 1 ' 1 ) , 
as obtained from a fit to the curve in Fig. 1. Figure 2 illustrates the main effects of the increase of 02 on RFS, showing 
that fcx5aFS f° r maximum growth is almost independent of 02, and the decrease in the range of unstable fc^. The 
increase of the growth rate with a 1: as predicted by Eq. (fT2")> . is also clear when we compare the growth rate for the 
monochromatic case, also shown in Fig. 2 (pw - plane wave), with the other scenarios. This scenario is not observed 
for RBS (cf. Figure 3); RBS shows not only the decrease of maximum growth rate with 02 but also the increase of 
the range of unstable wavenumbers, and the shift of the most unstable wavenumber as 02 increases, in very close 
agreement with kjf RBS . Furthermore, the lower bound of the range of unstable wavenumbers remains unchanged as 
02 increases since this limit is determined by fc™RBS' a function of o\. The magnitude of the growth rate is within 
the same order of magnitude for the full range of unstable wavenumbers which indicates, that depending of the noise 
source, the instability can easily grow in a wide range of wavenumbers. The trend of saturation of the growth rate 
with 02, as predicted for a waterbag distribution function, is also clear. 

We have studied electronic parametric instabilities in the presence of broadband fields; as expected, the maximum 
growth rates for these instabilities decrease with increasing 2 . Nevertheless, a qualitative difference has been identified 
between RFS and RBS in the underdense limit, and for 2 < fco, with RFS decreasing slower than RBS, thus 
less sensitive to broadband fields, illustrating a fundamental difference between these two processes. Numerical 
solutions of the generalized dispersion confirm the theoretical results and illustrate the increase of the range of 
unstable wavenumbers with 0. Generalization of this work to address the role of broadband radiation in self-focusing, 
stimulated Brillouin scattering is straightforward, thus allowing comparison with recent experimental results [l8| , and 
numerical particle-in-cell simulations. These results will be presented in future publications. 
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FIG. 1: Maximum growth rate as a function of the photon distribution width (02), normalized to the growth rate for 02 = 0. 
The photon distribution is defined by fco = 80ui p o/c, ao = 0.1 and o\ = 0.2 ko- 




0.997 0.998 0.999 

k L (U) p0 /C) 



FIG. 2: Growth rate of the instability in the range of &l leading to RFS, for a photon distribution with kh ko = 80ui p o/c, 
ao — 0.1 and tri = 0.2 fco, a - (72 = 0.29 fco, b - 02 = 0.39 fco, c - 02 — 0.59 fco, d - 02 — 0.78 fco, e - 02 = 0.98 fco- The plane 
wave (pw) pump scenario is obtained for the same parameters fco = 80lu p o/c, ao = 0.1 - RMI is also observed in this case. 
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FIG. 3: Growth rate of the instability in the range of kL leading to RBS for a photon distribution with ko = 80lo p q/c, ao = 0.1 
and <7i = 0.2&Q, a - <J2 = 0, b - <J2 = 0.25 ko, c - <72 = 0.5 ko, d - <72 = 0.75 ko, e - — ko- The range of unstable wavenumbers 
for the monochromatic scenario (pw) is also shown for reference, with a TaBSpw = 0.44 lu p q. 



